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SUMMARY 

V 

The  Program  Evaluation  and  Review  Technique  is  a 
very  well-known  and  widely-applied  method  for  estimating 
the  amount  of  time  required  to  complete  a  complex  project. 
However,  if  the  time  required  to  complete  each  subtask  is 
a  random  variable,  the  computation  of  the  expected  time 
to  complete  the  total  project  may  be  infeasible  (for  reasons 
of  dimensionality) .  The  simplest  way  to  estimate  the  ex¬ 
pected  time  to  complete  the  project  is  to  assume  each 
subtask  takes  exactly  its  expected  time  to  complete;  and 
then  find  a  critical  path.  This  procedure,  except  in 
very  simple  rases,  gives  an  underestimate  of  the  true 
expected  time  required.  In  1962  D.  R.  Fulkerson  published 
a  technique  for  estimating  the  expected  time  required 
which  gives  values  intermediate  between  the  simplest 
extimate  and  the  true  value.  This  Paper  suggests  another 
technique  which,  while  more  complex  than  Fulkerson's 
method,  is  more  accurate  in  that  it  gives  estimates  which 
fall  between  Fulkerson's  estimates  and  the  true  expected 
value. 


ANOTHER  ESTIMATE  OF  EXPECTED  CRITICAL  PATH  LENGTH 

IN  PERT  NETWORKS 

As  in  [1],  we  assume  that  we  have  a  PERT  network  with 
nonnegative  work  times  assigned  along  the  various  arcs  de¬ 
pending  on  a  discrete  random  variable  u,  where  th  probability 
of  the  value  u  occurring  is  Number  the  nodes  so  that 

j  <  i  if  there  is  an  arc  from  node  j  to  node  i.  We  also 
assume  that  work  times  on  arcs  coming  into  any  node  i  are 
independent  from  work  times  on  arcs  coming  into  a  different 
node  j.  In  other  words,  u  may  be  considered  to  be  a  Cartesian 
product  of  random  variables  u^  on  which  work  times  on  arcs 
coming  into  the  node  i  depend,  with  4r  (u)  being  a  product  of 
functions  <|k(u^). 

We  wish  to  approximate  the  expected  longest  path.  The 
simplest  approach  is  to  simply  assign  each  arc  its  expected 
work  time  and  calculate  the  longest  path  PQ  and  its  work  time 
for  this  assignment.  As  the  expectation  of  the  sum  of  random 
variables  is  the  sum  of  the  expectations  of  the  random  vari¬ 
ables,  this  equals  the  expected  length  of  a  critical  (longest) 
path  [3]  if  and  on? y  if  for  no  value  of  u  is  some  path  longer 
than  the  path  P  .  The  method  employed  in  (1]  has  greater 
flexibility  in  that  it  does  not  force  us  to  use  the  same 
path  for  all  events.  In  fact,  the  next  to  last  'ode  of  the 
path  is  allowed  to  depend  on  the  work  lengths  which  u  assigns 
to  arcs  to  the  final  node.  However,  for  each  u,  the  next  to 
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last  node  of  a  critical  path  also  depends  on  the  work 
lengths  which  u  assigns  to  earlier  arcs.  This  dependence 
is  the  reason  that  Fulkerson's  estimate  [1]  may  fall  short 
of  the  actual  expected  length  of  a  critical  path.  The 
following  example  illustrates  this. 


The  arcs  from  1  are  allowed  to  take  on  the  values  0,  1  inde¬ 
pendently  with  probability  £.  As  there  is  only  one  assign¬ 
ment  of  work  lengths  along  arcs  to  node  4,  Fulkerson's  method 
forces  us  to  go  through  node  2  all  the  time,  or  through  node 
3  all  the  time.  This  results  in  the  estimate  £  while  the 
actual  expected  critical  length  is  f.  The  problem  occurs  when 
for  a  fixed  assignment  of  work  times  on  arcs  to  the  final 
node  and  variable  assignment  of  work  times  to  earlier  arcs, 
there  is  a  close  decision  as  to  which  node,  j  or  k,  should 
be  the  next  to  last  node.  In  general,  andom  fluctuation 
will  make  a  path  through  j  longer  than  a  path  through  k 
for  some  u  and  vice  versa  for  other  u.  However,  if  a 
critical  path  to  j  tends  to  fluctuate  the  same  way  as  a 
critical  path  to  k  (a  critical  path  to  j  is  long  when  a 
critical  path  to  k  is  long)  then  still,  little  is  lost  in 
always  choosing  a  path  through  u  for  a  fixed  assignment  of 
work  times  along  arcs  to  i .  In  summary,  the  Fulkerson 
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estimate  is  weakest  when  there  is  a  large  variance  in 
critical  path  lengths  to  these  nodes  i  and  k,  and  there 
is  little  or  negative  correlation  between  these  critical 
path  lengths.  It  is  in  precisely  this  situation  that 
the  following  modification  aims  to  improve  Fulkerson's 
estimate. 


1  to  node  j.  We  define  the  following  quantities  by  recursion 

on  i  or  max  (i,j):  k^,  tu,  n^,  N^,  hLj  for  i  r  j,  and 

C. .  for  i  >  i.  Define  k,  =  h,  *  n,  =  N,  =  C. .  =  M.  .  =0. 

ij  J  1  i  1  1  ij  JJ 

Fix  i  and  let  S  be  the  set  of  nodes  from  which  there  are 

arcs  to  i.  Let  work  times  along  arcs  to  i  depend  on  the 

random  variable  t  and  work  times  along  arcs  to  nodes  j, 
j  <  i,  depend  on  the  random  variable  s  where  s  and  t  are 
independent,  the  probability  of  s  occurring  being  cp(s),  and 
the  probability  of  t  being  P( t) .  If  S  has  one  element,  let 
a(t)  =  b(t)  be  this  element.  Otherwise,  define  a(t),  b(t)  e  S 
so  that  for  j  e  S  -  {a(t),  b  ( t) } ,  hfl  ^  +  ta^(t)  -  hb(t) 

+  ^(t)^  —  ^  +  tj^t^  w*iere  tj(t)  is  the  work  time  along  the 
arc  from  j  to  i  for  the  random  variable  t.  In  the  following,  if 
a  term  has  0  in  its  denominator,  set  the  term  equal  to  0.  Let 

ki  "  l  P<C>  (ha(t)  +  ta(t)<t>)'  hi  “  rtP(t)  ^  (ha(t)  + 

ca  (t) (t) '  ^ha(t)  +  ta(t)(t)  +  hb(t)  +  tb(t)<t>)  +  K(L)]  whcre 

K(t)  =  [Cb(t)b(t)  ~  2Ca (t)b(t)  +  Ca (t)a (t)  +  (ha(t)  + 

Ca(t)(t>  ~  hb(t)  ~  cb(t)<t>>2I/2  max  (Mb(t)a(t)  + 

Lb(t)^  “  ta  ( t)  ^  ’  Ma(t)b(t)  +  ta(t)^t^  ~  ^(t)^^’ 
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ni  *  mtin  <na(t)  +  ta(c)<t»'  Ni  *  TX  [,"ax(Na(c)  +  ca(t)(c)’ 

^(t)  +  tb(t)(tW1,  for  1  >  ^Mji  “  <f*  <Mja(t)  " 

and  Mij  =  max  [max  (Ma(t)j  +  ta(t)(t),  +  tb(t;)  (t))  [, 

Cii  =  fp(t)  !Ca(t)a(c)  +  (ha(t)  +  ca(t)(t:)  ~  hi>2'  ~  2(hi  "  "i1 
(hi  -  kt),  and  for  i  >  j_,  =  Ep(t)Cfl -  k^  (Nj  -  tu). 

Define  C..  *  C. .  for  i  <  i.  We  shall  show  by  induction  on 
Ji  ij  J  7 

i,  that  for  each  node  i ,  there  is  a  function  y^u)  of  the 
random  variable  u,  satisfying  the  following: 

I  y^u)  <  xi»>  \  <  <  Ni,  Yj  (u)  -  yk(u) 

ya(t)(s)  +  ^(t)^  -  yi  (s>  ^  =  yi(u)  -  max  (ya(t)^s) 

+  Vc)(s)  +  cb(t)(c>); 

II  h.  -  y.  -  j:  Uu)y.  (u) ; 

u  x 

III  Ef(u) (y. (u)  -  hi)2  >  C. . ;  and 

u  L  1"L 

IV  ^(u)(yj(u)  -  hj)(yk(u)  -  hfc)  <  Cjk. 


Assume  that  the  above  statements  are  true  when  i,  j,  and 
k  are  replaced  by  integers  less  than  i.  For  each  s,  t  the 
longest  path  f-'-om  node  1  to  node  i  has  length  xj  t^(s)  + 
tj  (g  t)(c)  f°r  some  node  j(s,  t)  with  an  arc  from  node 
j(s,  t)  to  node  i.  Then  xj(s>  t)(s)  +  tj(s_  t)(c)  ;  xa(t)(s) 
+  ca(c)(t)-  Similarly,  x1(s>  t)(t)  +  tJ(s>  t)(t)  1  max 

(yd(c) <s>  +  ca(t)<c>-  yb(t)<s)  +  -  Hya(c)<s)  + 


5- 


ta(t)<t>  +  yb(t)<s>  +  tb(t)(t))  +^yb(t)(s)  +  cb<t)<c)  _ 
ya(t)<s)  -  ca(t)(t)  I  •  Then  ha  (t)  +  ‘atdW  ’  (ya(t)  (s) 

+  ca(t)(t))  -  §cp(s)max  lya(c)(s)  +  ta(t)(t),  yb(c)(s)  +  tb(t)(t)I 

*  i  |'i)(sHya(c)(s)  +  ta(t)(t)  +  yb(c)(s>  +  cb(t)(c)1  + 

i  ^9(s)|ya(t)(s)  +  ta(t)(t)  -  yb(t)(8)  -  cb(c)(c)| 

“  ilha(t)  +  Ca(t)(V  +  hb(t)  +  cbCc)  (t>  1  + 
i  ^(s)lyf(t)(s)  +  ta(t)(t)  -  yb(c)(s)  -  tb(t;)(c)|. 

Also,  i  rcp(s)|yb(c)(si  +  tb(t)(t)  -  ya(t)(s)  -  ea(t)(t>| 

>  S<p(8)[yb(t)(8)  +  tb(t)(t)  -  ya(t)(8)  -  ca(t)(t!]2/2  «x 

IMb(t)a(c)  +  Cb(t)(C)  ~  Ca(c)(C)’  Ma(C)b(t)  +  Ca(c)(t)  “  Cb(t)(C)' 
and  rM'Hyb(t)(s)  +  tb(c)(t)  -  ya(t)is)  -  ta(c)(t)]2 

=  E<f(si{  !(yb(C)(s)  -  hb(c))  -  (ya(t)(8)  -  ha(t))1  +  [(hb(t)  + 
Cb(t)(t>)  ~  (ha(t)  +  Ca(t)(r''12 

■  ^(s)(yb(t)(s)  -  hb(t))2  -  2  J*<»>(yb(C)<s)  -  hb(c)) 

(ya(t)(s)  "  ‘Vo’  +  M(sKya(c)(s)  -  ha(t))2  +  (hb(t)  +  tb(t)(t) 

“  (ha(c)  +  Ca(t)(t))'2 

-  °b(c)b(t)  “  2Ca(t)b(t)  +  Ca(t)a(t)  4  (hb(t)  +  cb(t)  - 
<ha(t)  +  ta(t)(t)l2' 

This  shows  that  we  may  define  yi(s,  t)  to  satisfy  the  following: 

ya(c)<s>  +  Ca(t)(t>  S  yi(s,  t)  S  max  <ya<t)(s)  +  Ca(t;(t)' 

yb(t)<s>  +  Cb(t)(t))  i  max  Ixa(t)<s)  +  xb(t)(s>  + 

tb(C)(t)]  <  xi(s,  t),  and  h.  -  y±  (t)v  (s)yi(s ,  t) 

1  ~  P(t)"  (s)x .  (s ,  t)  *  x,.  This  shows  by  induction  that  h.  lies 
s,t  i  i  1 

between  the  actual  expected  length  of  a  critical  path  and 
Fulkerson's  estimate  for  it. 
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Furthermore,  n£  =  min  [na(t)  +  ta(t)(C)]  <  min  Iya(t)<s> 

+  -  5J5  yi(s’  t)  -  j“*  yi<s,  t)  <  max  (max(ya^t)  (s)  + 

ta(t)^’  yb(t)^s^  +  tb(t)^t^ 

<  m^x  [max  (Ha{t)+  £a(t)(t),  Nb<t)  +  tb(c)(t»]  -  N.. 

For  i  >  j,  Mjl  <=  max  (Mja(t)  -  ta(t)(t» 

-  “?  lyj(s)  "  ya(t)<s>  "  ca(t)<c>>  2  iyj(s)  -  yi(£-  t)] 
and  Mi,  -  max  [max  <Ha(t)j  +  ta(t)(t),  ^(C)J  +  ^^(t))] 

>  max  [max  (ya(t)<«)  -  Yj(a)  +  ta(t)(t),  yb(t)(s)  -  y^a)  + 

"  “J  tmax<ya(t)(s)  +  ca(t)(t)’  yb(t)<s)  +  -  yj(s>] 

>  max  [  yi(s,  t)  -  (s) ). 

Also,  Cu  -  F.  P(t)ICa(t)a(t)  +  (ha(t)  +  ta(t)(t)  -  h.)2]  - 
2(hi  -  -  kt) 

<sFtP(t)m(S)[ya(|.)(s)  -  ha(t)  +  (ha(£)  +  ta(t)(t)  -  h.)]2 

-  2(ht  -  n^lhj.  -EtP(t)m(a)(ya(t)(s)  +  ta(t)(t))J 
-sitP(t)ct(sHya(t)(s)  +  ta(t)(t)  -  h^2 
-sFtP(t)m(s)(yl(8,  t)  -  ya(t)(s)  -  ta(t) (t)) (2hj  -  2nt) 
<sFtP(t)m(s)(ya{t.)(a)  +  ta(t)(t)  -  hl)2 
-s!:tP(tlm(.)(yl(s,  t)  -  ya(t)(s)  -  ta(t)(t))(2h1  -  y^s,  t) 

"  ya(t)<s>  -  ta(t)(t)) 

■  y.  P (t)m(s)  (y.  (s ,  t)  -  h.)2.  Finally,  for  i  >  i,  C.. 
s,t  1  1  1J 

■?  P(t)C.(t)j  +  <hi  "hj) 

2  r  fP(t)  ■  x(s)(ya(t)(s)  -  ha(c))(yj(a)  -  h.)J 

+sFtP(t)m(s)(yl(s,  t)  -  ya(t)(s)  -  ta(t)(t))(Nj  -  h,) 

>  i  iP(t)  ~  =p(s)(ya(t)(s)  -  ha(t))(y.(s)  -  hj)] 
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+gi:tP(t)cp(s)(yi(s,  t)  -  ya(t)(s)  -  ca(t) (c)) (y^ (s)  ~  hj) 

+  |  p (t)(ha(t)  +  ca(t)(t)  -  ht)  §  v(s)(y5(3)  -  hj) 

=sStp(t)cp(s)(yi(s,  t)  -  hi)(y^(s)  -  h^). 


h4=i(i+0+^+0)+(A-0+i  +  (^+0-^  -  0)2)/2  -  3/4, 
the  actual  expected  critical  length.  Naturally,  the  example 
exaggerates  the  accuracy  of  our  estimate.  One  problem  is 
that  for  some  assignment  of  work  times  to  arcs  into  a  node 
i,  there  may  be  a  close  decision  between  three  or  more  nodes 
as  to  which  should  iranediately  precede  i  on  a  critical  path. 

In  the  example 


our  estimate  is  still  3/4,  but  the  actual  expected  critical 
path  length  is  1  -  2  n.  Without  using  the  assumption  that 
work  times  on  arcs  in  bundles  leading  into  nodes  j,  j  <  i, 
are  independent,  it  is  not  possible  to  use  the  correlations 
between  critical  path  lengths  coming  in  two  nodes  to  calcu¬ 
late  how  much  critical  paths  through  c(t)  to  i  increase 
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™ax  (ya(t)(a)  +  ta(c)(t),  yb(t)(s)  +  tfc(t)(t);  that  is, 
by  how  much  g  cC(s)max  (ya(t)(s)  +  ta(t)(t)y  yb(t)(s)  + 

^(t)^*  yc(t)(s)  +  exceeds  |  cp (s)max  (ya(t)(s)  + 

ta(C)^t^  yb(t)^s^  +  Vt)(t)>*  In  fact>  suppose  that 
ta(l)(l)  =  Cb(l)^  =  Cc(l)^  =  ®  and  s  ta^es  on  f°ur  values 
with  probability 


xa(l)<s)  ’ 

xb(i)(s)  - 

xc(l)<s>  * 


s  : 

?a(l)(s> 

Vl)(s) 


'c(l) 


(a)  - 


I 

0 

0 

0 


II 

0 

1 

1 


III 

1 

0 

1 


IV 

1 

1 

0 


Then  there  is  no  correlation  between  any  two  of  the  follow¬ 
ing:  ya(t)(s)>  yb(i)(s>'  and  yc(i)(s)*  However,  it  is  never 
necessary  to  use  c(l)  in  a  critical  path;  regardless  of  s, 
there  is  a  critical  path  through  a  or  b.  However,  in  some 
sense  there  is  correlation,  as  /nN  =* 

a  (l)b(l)c (1) 

r  :p(s) (ya(i) (s)  -ya(i))^b(i)<8>  -  ^(D^ycd)^'  -  yC(i)) 

-  H  -  i>(  -  i)(  -  £)  +  3(J)(i)(4)(  -  -  1/8.  There  are 

two  problems  in  trying  to  use  this.  First,  we  would  have  to 

calculate  an  expression  for  each  triple  of  nodes. 

Secondly,  we  would  have  the  problem  of  using  a  third  degree 

term  C.  with  the  second  degree  terms  C... 

We  have  already  experienced  these  difficulties  to  some 

extent.  We  have  to  calculate  expressions  C. .  and  M. .  for 

ij  iJ 

each  pair  of  nodes,  whereas,  in  Fulkerson's  estimate,  it 
was  only  necessary  to  calculate  for  each  node.  Another 
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problem  was  that  we  were  able  to  calculate 

§  "(s)(yb(c)(s)  +  T,(t)(t)  ~  >'a(t)<s>  ~  but  we 

needed  |  „(s)  |yb(c)  (■)  +  t^^t)  -  ya(c)(s)  -  ta(t)(t)| 
instead.  Our  solution  was  to  bound  the  ratio  by  the 

■*>jlyb(t;)(s>  +  tb(c)(t)  -  ya(t)(s)  -  ta(t)(t)|.  With  tiore 

simplicity  and  less  sharpness,  we  could  have  used 

max  (Nb(t)  +  cb(c)  ~  na(t)  “  Na(t)  +  ta(t)^t) 

~  "b(t)  "  cb(t)(t»  instead  °£  Mx  (Mb(t)a(t)  +  tb(t)<t) 

”  ca(t)<c)-  Ma(t)b(t)  +  ~  ^(t)^)'  as  a  bound  £or 

this  ratio.  Either  of  these  bounds  does  not  lose  very  much 

sharpness  if  yb(t)(s)  +  tb(t)(t)  -  ya(t)(s)  -  tfl(t)(t)  as 

a  function  of  the  random  variable  s  has  a  high  probability 

of  taking  on  a  value  near  the  endpoints  of  its  range,  as  in 

our  example.  However,  this  will  not  usually  be  the  case, 

particularly  if  we  have  a  continuous  random  variable  or 

if  the  node  i  is  far  from  node  1  so  that  lots  of  random 

fluctuations  make  yb^(s)  +  tb(t)(t)  ~  ya(t)(s)  “  ta(t)^t) 

behave  nearly  like  a  continuous  distribution.  One  solution 

is  to  assume  that  yb(t)(s)  +  tb(t)(t)  -  ya(t)(s)  -  ta(t)(t) 

is  distributed  for  fixed  t  as  a  normal  distribution.  If 

the  random  fluctuations  add,  then  this  is  a  very  plausible 

assumption,  if  the  node  i  is  far  from  node  1,  by  the  central- 

limit  theorem.  Now  yb(t)(s)  +  tb(t)(t)  “  ya(t)(s)  “ 

has  mean  x  «  hb(t)  +  tb(t)(t)  -  ha(fc)(t),  -  tfl(t)(t)  and 
2 

a  -  Cb(t)b(t)  “  2Ca(t)b(t)  +  Ca (t)a (t) *  Let  r  “ 

ha(t)  +  ta(t)(t)  "  hb(t)  ~  ^(t)^^  Then  e  cp(s)  |yb(t)  (s)  + 

^(t)^  “  >'a(t)<s)  “  ta(t)<t>  I 
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approximately  equals 


CO 

Jj*b(t)(s)  +  tb(t)<c)  ~  VO(s)  -  ca(t>(c>  I  yfea 


r  !yb(c)(s>  -  hb(t)  -  ya(t)<s)  +ha(t)>2/2°2 


d<*b(t)<s)  -  ya(t)<s» 


-E 


1  (x  -  r)  a"  *2/2“2dx 


r  J7na 


+  fr  -=r  (r  -  x)  e  x  dx 
J-to  v^tto 


2  a 


r2/2a2  +  2r _  rr/a  g  -  xz/2 

7*T  Jo 


ysr 


dx 


where  the  last  term  can  be  expressed  in  terms  of  the  error 
function: 


i  pr/a  -  x2/2 

JZ?  J0  e  dlt  ■  4  ~  err(r/o) 


The  expression 
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f(r,a)  *  ~pF~~  e 

«/2rr 


2, ,  2 


lla  r  ™  l*r/a  e  x  /2dx  if  a2  >  0, 
/2rr  °  o 


_•*  2 


f(r,cr)  =  y  >0  if  oZ<  0,  where  r  -  ha(>^  +  ~ 


a(t)  T  a  (t) 


hb(t)  tb(t)(t)  and  a  "  Cb(t)b(t)  2Ca(t)b(t)  +  Ca(t)a(t) 


may  be  used  to  replace  K(t)  in  recursively  defining  h.^.  Note 
that 

M_  ,  _i _  e  -  r2/2a2  .  _£ _  e  -  r2/ 2a2  _  r2  -  r2/2a2 

30  2JT5t  a  2J& 

,  _L_  e  -  r2/-"2  >  0. 

75T 

Therefore,  f(r,  a)  increases  as  a  increases,  and  f(r,  a)  > 

lim  f(r,  cj)  =*  y-  *=  f(r,  0).  Therefore,  we  may  define 
a-*0  z 

hi  *  i  p(t)^ha(t)  +  ta (t)  +  hb(t)  +  ^(t)^^  +  £^r>  a)^' 


Also,  our  conservative  estimates  of  a  give  a  conservative 
estimate  for  f(r,  a).  Advantages  of  the  use  of  f(r,  a)  are 
closer  estimates  of  the  expected  critical  patn  lengths  and 
avoidance  of  the  use  of  the  expressions  and  M^.  The 
major  disadvantage  is  that  we  no  longer  have  a  guarantee 
that  our  estimate  for  the  expected  critical  path  length,  will 
be  optimistic,  although  it  usually  will  be.  In  fact,  consider 
the  example 
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where  the  arc  from  1  to  3  takes  on  the  work  time  0,  p  of  the 
time;  1,  p  of  the  time;  and  |  the  rest  of  the  time.  Then 
^22  =  C33  ~  i(2p)*  The  actual  expected  critical  path 

length  is  ^  +  p/2.  The  new  estimate  for  h^  is  |(|  +  0  +  |  +0)  + 


As  p  gets  small,  sooner  or  later  this  exceeds  the 


expected  critical  path  length. 


Even  when  there  are  too  many  nodes  to  keep  track  of  the 


for  each  pair  of  nodes,  it  is  still  possible  to  carry 
through  this  procedure  for  the  first  hundred  or  so  nodes. 


label  the  bu  as  f^,  and  continue  with  Fulkerson's  recursive 
definition  [1]  of  the  f^.  When  our  original  method  of  defin¬ 
ing  the  h^  is  used,  we  still  obtain  an  estimate  which  is  no 
larger  than  the  expected  critical  path  length. 
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